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ABSTRACT 

A generalized Newtonian potential is derived from the geodesic motion of test particles in 
Schwarzschild spacetime. This potential reproduces several relativistic features with higher 
accuracy than other commonly used pseudo-Newtonian approaches. The new potential repro- 
duces the exact location of the marginally stable, marginally bound, and photon circular orbits, 
as well as the exact radial dependence of the binding energy and the angular momentum of 
these orbits. Moreover, it reproduces the orbital and epicyclic angular frequencies to better 
than 6%. In addition, the spatial projection of general trajectories coincide with their relativ- 
istic counterparts, while the time evolution of parabolic -like trajectories and the pericentre 
advance of elliptical-like trajectories are both reproduced exactly. We apply this approach to a 
standard thin accretion disc and find that the efficiency of energy extraction agrees within 3% 
with the exact relativistic value, while the energy flux per unit area as a function of radius is re- 
produced everywhere to better than 7%. As a further astrophysical application we implement 
the new approach into a smoothed particle hydrodynamics code and study the tidal disruption 
of a main sequence star by a supermassive black hole. The obtained results are in very good 
agreement with previous relativistic simulations of tidal disruptions in Schwarzschild space- 
time. The equations of motion derived from this potential can be easily implemented into 
existing Newtonian hydrodynamics codes at virtually zero additional computational effort. 



1 INTRODUCTION 



Our current understanding of some of the most energetic phenom- 
ena in the Universe (such as active galactic nuclei, X-ray binar- 
ies and gamma-ray bursts) involves the accretion of gas onto as- 
trophysical black holes as the underlying mechanism for power- 
ing these sources (see e.g. prank et al.||2002| >. It is clear that the 
proper study of any of these systems should account in a consistent 
way for the strong gravitational field that prevails in the vicinity of 
these relativistic objects. It is nonetheless remarkable that many of 
the early works on the subject, which were essentially Newtonian 
with few general relativistic effects incorporated 'by hand', proved 
to be very successful at modelling a variety of accreting systems. 
For instance, the standard thin disc model of Shakura & Sunyaev] 
{1973} is purely Newtonian and the only result from general re- 
lativity that it uses is the existence of the marginally stable circular 
orbit at which their disc model was truncated. Similarly, the model 
of a thick accretion disc introduced by Paczyhsky & Wiita {1980} 
(PW hereafter) is based on Newtonian dynamics with the substitu- 
tion of the Newtonian gravitational potential $ N = —GM/r with 
the pseudo-Newtonian potentiaQ 

$pw(r) = GK 0-, (1.1) 

r - 2r g 

where r g = GM/c 2 is the so-called gravitational radius. The po- 
tential $ PW not only reproduces the correct location of the margin- 

1 The qualifier 'pseudo' is used to indicate that the related potential does 
not satisfy the Poisson equation V 2 <1? = AirGp. 



ally stable and marginally bound circular orbits around a Schwarz- 
schild black hole (at r m5 = 6r g and r mb = 4r g , respectively) 
but it also gives reasonable approximations to other quantities, 
e.g. binding energy (percent error (p.e.) 13%), orbital frequency 
Q, (p.e. ^ 50%), epicyclic frequency fl" (p.e. ^ 84%), and peri- 
centre advance (p.e. ~ 30%). This potential has been used in a 
large number of studies of accretion flows onto non-rotating black 
holes to mimic essential general relativistic effects within a New- 
tonian framework (e.g. Matsumoto et al. 1984, Abramowicz et al. 
T988l|Chakrabarti & Titarchuk|1995||MacFadyen & Woosley|1999| 
Hawley & Balbus 2002 , Lee & Ramirez-Ruiz 2006, Rosswog et al. 



2009| >. 



Other pseudo-Newtonian potentials have been introduced to 
give better approximations to specific general relativistic features 
but at the price of reproducing some other properties with less ac- 
curacy. For instance, the pseudo-Newtonian potential introduced by 
|Nowak & Wagoner|(l99T] > (NW hereafter) 

GM / 3r g . 12 



$Nw(r) = 



+ 



(1.2) 



reproduces the angular frequencies SI and fi" with better accur- 
acy (p.e. ^ 14% and ^ 42%, respectively) than $ PW while 
still giving the correct location of r m3 . However, it locates r mb at 
~ 3.5 r g and gives a less accurate estimate of the pericentre ad- 
vance (p.e. > 80%). See |Artemova et al.|(l99oT > for a comparison 
of the performance of $ PW , $»w, and other pseudo-Newtonian po- 
tentials in reproducing the structure of a thin disc around rotating 
and non-rotating black holes. 

The fact that no pseudo-Newtonian potential can repro- 
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duce equally well different dynamical properties of Schwarzschild 
spacetime has limited their range of applicability. For instance, a 
poor estimation of the orbital and epicyclic frequencies will hamper 
the ability of capturing time-dependent behaviour, such as the onset 
of instabilities in an accretion disc that might eventually lead to ob- 
served signatures (such as quasi periodic oscillations, see e.g. |Kato| 
|2001[ >. On the other hand, a large error in the estimation of the bind- 
ing energy of Keplerian circular orbits will lead to inaccurate estim- 
ates of the total luminosity of accretion discs. Another issue that 
is usually forsaken is that most pseudo-Newtonian potentials are 
designed for accurately reproducing circular orbits but not neces- 
sarily more general trajectories. Nevertheless, they are frequently 
used in applications in which correctly reproducing general traject- 
ories might be of crucial importance (e.g. the collapsing interior of 
a massive star towards a newborn black hole or successive passages 
of a star orbiting a black hole before becoming tidally disrupted). 

In the present work we propose a generalization of the Newto- 
nian potential that accurately describes the motion of test particles 
in Schwarzschild spacetime. Even though our approach is formu- 
lated in an entirely Newtonian language, it is, indeed, equivalent 
to the one of Abramo wicz et all] {1997), in which the authors de- 
scribed particle dynamics under Newtonian gravity in a curved 
space. The main difference is that, instead of modifying the back- 
ground geometry, we modify the gravitational potential. In addition 
to the Newtonian 1 /r-term, our potential includes an explicit de- 
pendence on the velocity of the test particle. Following the nomen- 
clature of Lagrangian mechanics (see e.g. |Goldstein et al~| 2002), 
we therefore call it a generalized Newtonian potential [^Moreover, 
the present potential is derived from the actual geodesic motion 
of test particles in Schwarzschild spacetime. This is in contrast to 
most pseudo-Newtonian potentials that are introduced as ad-hoc re- 
cipes or as fitting formulae that mimic certain relativistic features 
(see, nevertheless Abramowicz 2009 for a a posteriori derivation 

of $pw). 

The remainder of the paper is organised as follows. In Sec- 
tion[2]the generalized potential and the corresponding equations of 
motion are derived and contrasted with the exact relativistic expres- 
sions. In Section[3]we compare the performance of our generalized 
potential with $ N , $ PW , and $ NW in reproducing several dynam- 
ical features of test particle motion in Schwarzschild spacetime 
such as purely radial infall, Keplerian and non-Keplerian circular 
motion, general trajectories, pericentre advance, and two simple 
analytic accretion models (the thin disc model of Shakura & Sun 



yaev|1973|as generalized to Schwarzschild spacetime by 



Novikov 



& Thorne|1973| and the toy accretion model of Tejed a et al.|2012| )" 
At the end of this section we implement the new potential into a 
Newtonian smoothed particle hydrodynamics (SPH) code and sim- 
ulate the tidal disruption of a solar-type star by a supermassive 
black hole. Finally, we summarise our results in Section]?] 



2 GENERALIZED NEWTONIAN POTENTIAL 

Consider a test particle with four-velocit)|^]ii f ' = dx^/dr follow- 
ing a timelike geodesic in Schwarzschild spacetime (where r is 
the proper time as measured by a co-moving observer). Given the 
staticity and spherical symmetry of the Schwarzschild metric, the 
motion of the particle is restricted to a single plane (orbital plane) 
and it is characterized by the existence of two first integrals of mo- 
tion: its specific energy £ and its specific angular momentum h s 
given by (e.g. |Frolov & Novikov|1998) 



£ = 



-Ut = c 



_ 2r, 
r 

2 dip 



frr' 



(2.1) 



(2.2) 



where ip is an angle measured within the orbital plane. These two 
equations together with the normalisation condition of the four ve- 



locity, u^u 11 
tion 

dr 

d7 



-c , lead to the equation governing the radial mo- 



£- 



+ 



2 GM hi 



2r, 



(2.3) 



Using Eq. {2A\ , we can rewrite Eqs. \2.2\ and < |2. 3| > in terms of de- 
rivatives with respect to the coordinate time t, i.e. 



dr 
dt 

where 



c 

J 



1 - 



dp 
~dt 

2r e 



c 

J 



2r, 



(2.4) 



2E S + 



2 GM h 2 



1 - 



2r 



E s = 



2 c 2 



(2.5) 



(2.6) 



This is a suitable definition of energy since, in the non-relativistic 
limit (nrl) in which r/r g 3> 1 and v 2 /c 2 <C 1 (where v can be 
either r or rip and the dot denotes differentiation with respect to the 
coordinate time i), it converges to the specific Newtonian mechan- 
ical energy, i.e. 



jp . 771 / -2 | 2 -2\ GM 

E s — > E N = - (r + r ip ) . 

nrl 2 v ' r 



(2.7) 



Our starting point for a generalization of the Newtonian po- 
tential is E s in the limit where v 2 /c 2 <C 1, but not necessarily 
r/r g ^> 1, i.e. 



E s 



-> E G 



(r-2r g 



+ 



2r E 



GM 
r 



This equation can be recast as 

E G = T + <f> G (r,f,ip), 



(2.8) 



(2.9) 



where T = (r 2 + r 2 ip 2 ) /2 is the non-relativistic kinetic energy 
per unit mass and $ G is a generalized Newtonian potential given 
by 



2 The potential introduced by Semerak & Karas ( 1999J for a Newtonian 
description of test particles motion around a rotating black hole is another 
example of such a generalized potential. However, when this potential is 
applied to the non-rotating case, it does not give an overall satisfactory per- 
formance (e.g. one finds r mb = r ms = 0). 



<f> G (r,-r,ip) 



GM 



2r„ 



r - 2r, 



r - 2r, 



.2 , r 2 p 2 
r + 



2 

(2.10) 



3 Greek indices run over spacetime components, Latin indices run only 
over spatial components and the Einstein summation convention over re- 
peated indices is adopted. 
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The first term on the right hand side of Eq. ( |2,10^ is the usual New- 
tonian potential that dictates the gravitational attraction due to the 
interaction of the central mass and the rest mass of a test particle, 
while the second term can be interpreted as an additional contribu- 
tion due to the kinetic energy that is also gravitationally attracted 
by the central mass. Contrary to a pseudo-Newtonian potential, $ G 
does satisfy the Poisson equation (when the only source of the grav- 
itational field is the central mass). 

We now use $ G to construct the following Lagrangian (per 
unit mass) 



L = T - $r 



2-2 3-2 

r r r ip 

+ ^~ 

r — 2 r„ 



2r g ) 2 



+ 



GM 



(2.11) 



Since L is independent of t, the energy E G as defined in Eq. l |2.8| > is 
indeed a conserved quantity of the corresponding evolution equa- 
tions. On the other hand, the independence of L on ip guarantees 
that the specific angular momentum defined as 



dL 



r 3 if 
-2r. 



(2.12) 
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Figure 1. Infall time for a test particle falling radially from infinity as calcu- 
lated in general relativity and using the potentials <E> N , "^g, 'fpw, and 3>nw- 
In all the cases we have taken = —0.3 c and have chosen r = 20 r s as 
synchronization radius. Note that both Ate and Ats diverge to infinity as 
the particle approaches the black hole horizon located at r = 2 r s . 



is also conserved. By combining Eqs. < |2. 8| > and (2.12| l we get the 
following expression for the radial motion 



dr 
dt 



2r, 



2E G + 



2 GM 



It 2 

»G 



2r« 



(2.13) 

which has a clear resemblance to the exact relativistic expression 
in Eq. (23} and, as we shall show below, correctly reproduces a 
number of relativistic features. Note in particular that for particles 
for which £ = c 2 (i.e. parabolic-like energies), Eqs. ( |2.12} and 
l |2.13F are identical to their relativistic counterparts (Eqs. |2.4| and 
2.5| respectively) once the correspondences 



E G ^ E s 



and 



h G -H> h s 



(2.14) 



have been taken. 

Even though the whole evolution of the test particle motion 
is already determined by Eqs. l |2.12[ l and l |2.13fr , it is also useful to 
compute the corresponding expressions for the accelerations com- 
ing from the Euler-Lagrange equations for an arbitrary coordinate 
system (r, 9, </>), i.erl 



GM 

r2 



1 - 



2r E 



+ 



v \ r J r(r — 2 r g ) 
+ (r-3r g ) (e 2 + sin 2 ° " A 
2r0 (r- 3r E 



r 
r 



r - 2r g 
r — 3 r. 



2r B 



+ sm cos f 
— 2 cot 8 d> 



(2.15) 
(2.16) 
(2.17) 



These equations should be compared against the general relativistic 



The angular velocities 9 and ej> are simply related to ip by the relation 



ones which are given by 



GM 



2r« 



2r E f 2 



2 -4 
£ 2 r(r — 2 r g ) 



) 2 + sin 2 l 



2f0 
r 

2r4> 
r 



r - 3r g 
r - 2r g 
r — 3 r t 
r — 2 r. 



+ sin 6 cos f 



2 cot e < 



(2.18) 
(2.19) 

(2.20) 



from where we can see that they are identical except for the factor 
c 4 /£ 2 multiplying the first term in Eq. ( |2.18^ . 

For an implementation of the present approach into an ex- 
isting hydrodynamics code, the acceleration components may be 
needed in Cartesian coordinates. The corresponding expressions 
are provided in Appendix |A| 



3 COMPARISON WITH PREVIOUS APPROACHES 

In the following subsections we compare the performance of $pw, 
$ NW , and $ G (Eqs. |l.l| |1.2| and |2.10| respectively) in reprodu- 
cing several relativistic features of the motion of test particles in 
Schwarzschild spacetime. As a reference to illustrate the import- 
ance of relativistic effects, we also show the results obtained by 
applying the Newtonian potential <& N . 



3.1 Radial infall 

Consider a particle in radial free-fall, i.e. <j> — 9 — 0. From 
Eq. l|2.13|) it follows that the amount of time that it takes for the 



particle to fall from a radius r 2 to a smaller radius T\ is given by 

2r E \ [~ 2 GM 



At c 



2E G + 



dr. 



(3.1) 



It is clear that At G will coincide with the corresponding relativistic 
value At s (as calculated from Eq. |2.5fr only when E s = E G = 0, 
i.e. for a vanishing radial velocity at infinity. In the general case we 
will have At G > At s if E G < and At G < At s if E G > 0. Also 
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Figure 2. Comparison of different quantities associated to circular motion. From top to bottom, left to right the panels show: the specific energy E c , the specific 
angular momentum h", the orbital angular velocity Q, and the epicyclic frequency for small perturbations parallel the orbital plane fill . Analytic expressions 
for the different quantities plotted in this figure are collected in Appendix [B] 



note that both At G and At s diverge as r x — > 2 r s , which coincides 
with the description made by an observer situated at spatial infinity 
in Schwarzschild spacetime. In Fig. [T] we show the inf ailing time as 
calculated from the relativistic solution and compare it with those 
coming from the use of <3> G , <&pw, and $ NW - 



3.2 Circular orbits 

We consider now the special case of circular orbits as defined by 
the conditions 



r = and r = 0. 



(3.2) 



After substituting these two conditions into Eqs. q2TT3j and q2TT5j 
we get a system of two equations that can be solved for the corres- 
ponding values of h G and E G as 



(3.3) 

E%=-^(^A, (3.4) 



V r - 3r E ' 
GM (r -4r t 
2 r \ r — 3 r f 

which are identical to the relativistic expressions and lead to the 
exact locations for the corresponding photon orbit r ph = 3 r g (at 
which both h G and E G diverge), marginally bound orbit r mb = 4 r g 
(for which E G = 0), and marginally stable orbit r ms = 6r g (at 
which both h G and E G reach their minima). See the top panels of 



Figure [2] for a comparison of E c and h c as calculated using the 
different potentials. 

On the other hand, by combining Eqs. ( |2.12^ and ([33} one 
finds the orbital angular velocity 



/ GM 


fr-2r s \ 


1 r-3r s 


V ^ J 



(3.5) 



that should be compared against the exact expression in Schwarz- 
schild spacetime 



GM 



(3.6) 



For r 6r g , Q G reproduces the exact value with an accuracy better 
than 5.7%. In the bottom left panel of Figure [2] we compare these 
two frequencies as well as the ones corresponding to the potentials 
$n, $pw, and $nw- The analytic expressions of all of the different 
quantities plotted in this figure are summarised in Appendix [B| 



3.3 Perturbed circular orbits 

We now calculate the epicyclic frequencies associated with a small 
perturbation of a stable circular orbit in the equatorial plane. The 
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unperturbed trajectory satisfies 

x\t) = (r, tt/2, fJ G t), 
x\t) = (0,0, fi G ), 

while for the perturbed trajectory we have 

x'\t) = (r + 5r, tt/2 + 89, Q G t + 5<j>), 
x'\t) = (5f,5e,fl G +64>). 



(3.7) 



(3.8) 



Substituting these expressions into Eqs. ( |2.15F - {2.17| > we obtain the 
following system of (linearised) differential equations for the per- 
turbed quantities 



5r 



n , + 2GM (r _ 2rg)(r _ 



■ 4r„ 



5r 



+ 2fi G (r-3r g 



2Q C 
v 



3 r { 



Sr. 



(3.9) 
(3.10) 

(3.11) 



From Eq. j3.10| > we see that, to first order, the vertical perturbation 
decouples from the other two directions and has the same angular 
frequency as the orbital motion, i.e. 



Or- — £lc 



(3.12) 



as is also the case for Schwarzschild spacetime. On the other hand, 
from Eqs. \3.9\ and {344} we can se e that the radial and azimu thal 
perturbations are coupled. Following 



Semerak & Zacek 



(20001, we 



assume that the solution to both equations is a harmonic oscil- 
lator with common angular frequency $1 G , i.e. 5r — 5r e 1 °g * and 
5(f) — 8<j) e lSl G t where 5r and 5(f> are constant amplitudes. After 
substituting these solutions into Eqs. $.9\ and l |3.1 1\ we get as the 
only non-trivial solution 



GM 



6r, 



3r s 



that should be compared against the exact relativistic value 



GM 



6 r. 



(3.13) 



(3.14) 



Again, for r 6 r s , Q G reproduces the exact value with an accur- 
acy better than 5.7% (see bottom-right panel of Figure[2}. 

Note that, even though none of the three frequencies Q G , fl G 
and qJI coincides with the corresponding relativistic expressions, 
they do keep the same ratios, i.e. 



G 

fie 



(3.15) 



— ^ = 1, and — - = — - -■ 

\IS W G iiS 

The p seudo-Newtonian potential introduced by Kluz niak & Lee| 
(2002 i also reproduces exactly the ratio fl"/fl, although it does 
not reproduce satisfactory other important features (e.g. one ob- 
tains r ms = 3 r g , r mb = 1.3 r g , and a p.e. of ~ 42% for the binding 
energy). 

3.4 Non-Keplerian circular orbits 

Consider a test particle going round a circular orbit with uniform 
angular velocity = d(j>/dt. This velocity does not necessarily co- 
incide with the Keplerian value and, in order to sustain this circular 
trajectory, an external force must be continuously applied along the 
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Figure 3. Radial thrust needed to keep a rocket hovering at a fixed radius. 
Note that J-q qualitatively follows the relativistic value of f and vanishes at 
the event horizon (r = 2 r g ). This is not a physical result but rather a reflect 
of the fact that the Schwarzschild coordinates are ill-behaved at this radius. 



radial direction. The actual nature of this force (e.g. pressure gradi- 
ents, electromagnetic or propulsion from a rocket) is unimportant 
for the present discussion. We can calculate the necessary radial 



thrust from Eq. l |2. 15| > as 

GM 



1 - 



2r. 



— Q (r — 3 r g 



(3.16) 



On the other hand, the corresponding relativistic expression is given 
by 

-1/2 



Fs = a = 1 - 



2r- 



(3.17) 



where a r is the physical value of the radial component of the the 
four-acceleration a M = u"u f ' ; „ (i.e. as projected along the radial 
direction of a local tetrad). Writing the four-velocity of the test 
particle as = it*(l, 0, 0, fl), it is simple to check that the only 
non-vanishing component of the four-acceleration is 

2rA GM/r 2 — rQ 2 ^ 



l-2rJr-r 2 Q 2 /c 2 '' 
from where we can rewrite Eq. l |3.17| > as 



GM 



n 2 



3r„ 



2rjr -r 2 Q, 2 /c 2 



1 - 



2r« 



-1/2 



(3.19) 

In Figure [3] we compare T as obtained from the different po- 
tentials for f2 = 0, i.e. a rocket that hovers still at a fixed radius. 
From this figure we see that J- G gives the worst approximation to 
the relativistic value (p.e. ^ 64% for r > 6r g ). In particular note 
that J-s diverges to infinity as r — > 2 r g while _7- G vanishes at this 
radius. However, the behaviour of F G is in qualitative agreement 
with the relativistic expression forf (= GM(r - 2r g )/r 3 ) which 
is a clear indication that the description made using <!? G does not 
correspond to a local observer but rather to one located at spatial 
infinity. 

Note that Eq. {346} captures an important feature of the exact 
relativistic expression in Eq. {349}, namely that the radial thrust 
becomes independent of the angular velocity at the photon orbit 
r = 3r g ( [Abramowicz & Lasota||1986fr, and that the centrifu- 



gal force reverses sign at this radius l [Abramowicz| 1 990} . As | Ab 
|ramowicz~& Miller ( 19901 pointed out, a consequence of this ef- 
fect is the general relativistic result that the ellipticity of a slowly 
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Figure 4. Comparison of an elliptic-like trajectory in Schwarzschild spacetime with those obtained using <3> N , <J> G , ^pw. an d $kw- In all the cases the 
pericentre and apocentre of the orbit have been fixed as r p = 5 r s and r a = 40 r g , respectively. The common initial point is indicated by an open circle. The 
left panel shows the spatial projection of the different trajectories onto the x-y plane (with x = r cos tp and y = r sin tp). The panels on the right hand side 
show the corresponding azimuthal (top) and radial (bottom) velocities as function of the radius. Note that the spatial projection of the trajectory associated to 
the potential $ G reproduces exactly the relativistic one, while the velocities are very well approximated. 



rotating Maclaurin spheroid changes monotonicity at r « 5r g 
as the spheroid contracts with constant mass and angular mo- 
mentum (whereas in the Newtonian case the ellipticity monotonic- 
ally decreases as the mean radius of the shrinking body decreases) 



(Chandrasekhar & Miller 1974 


MillerM977). Following the same 


approach as Abramowicz & Mi 


ler 1 1990 1, it is simple to check that 



the present Newtonian description leads to the same expression for 
the ellipticity as their Eq. (12') that quantitatively reproduces this 
effect (p.e. < 35% for r > 4r g ). 



3.5 General orbits and pericentre advance 

The geometric description of the trajectory followed by a test 
particle is obtained by combining Eqs. ( |2.12} and < |2. 1 3^ as 



dtp 



he 



where 



P(r) = 2E G r i + 2GMr s - h 2 G r{r 



2r K 



(3.20) 



(3.21) 



Eq. ( |3,20F is formally identical to the general relativistic expres- 
sions once the correspondences in Eq. (2.14| l have been taken. This 
means that the spatial projection of the trajectories obtained with 
$ G is identical to that coming from the full Schwarzschild solution. 
In other words, the orbital parameters (e.g. apocentre, pericentre, 
eccentricity) and characteristics (e.g. pericentre advance, whether 



it is bound, unbound or eventually trapped) will have the same 
functional dependence on the corresponding constants of motion. 
In Figure|4]we show an example of a generic elliptic-like trajectory 
as resulting from the full general relativistic solution and compare 
it with those coming from the use of <I> N , $ G , <&pw, and $nw- 

Just as in the general relativistic case, the solution to Eq. < |3.20| > 
can be written in terms of elliptic integrals (see e.g. Chandrasekhar] 
|1983| [Tbjeda et aL||2012) . In the particular case of a bound orbit 
(elliptic-like trajectory) such that r £ [r p , rj (where r p and r a are 
the pericentre and apocentre of the orbit, respectively) we define its 
(half) orbital period as 



n c 



h G dr 



2h G K(k) 



\/PM ^-2E G r p (r a -r b )' 



(3.22) 



where K(k) is the complete elliptic integral of the first kind whose 
modulus k is given by 



r b (r a — r p ) 
r p (r a — r b ) ' 



(3.23) 



and r b = GM/E G — r a — r p is the smallest positive root of P(r) 
The pericentre advance or precession is given by 



pericentre advance = n G — it. 



(3.24) 



It is simple to check that, in the non-relativistic limit, r b — > 
and, thus, k — > 0. Given that K(0) = tt/2, it follows then from 
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Figure 5. Pericentre advance as a function of the eccentricity e for a fixed 
apocentre r a = 40 r s . Note that <J> G reproduces exactly the relativistic 
value, while the result obtained with <f> PW is consistently below the exact 
value by ~ 30%. On the other hand, for most values of e, "J>nw produces a 
negative shift. 



Eq. ([3T22J that 



n c 



(3.25) 



In Figure[5]we have plotted the pericentre advance as a func- 
tion of the eccentricity e = [r a — r p )/(r a + r p ) and as predicted by 
general relativity, "1> N , <J> G , $pw, and $ NW . As already mentioned, 
the pericentre advance resulting from $g corresponds to the exact 
relativistic value. On the other hand, the pericentre advance ob- 
tained from <5>pw is off by ~ 30%, while, for most values of the 
eccentricity, $ NW yields a negative shift (i.e. a pericentre lag rather 
than advance). 

The pseudo-Newtonian potential introduced by Wegg (2012) 
was specifically designed to give an accurate description of the 
pericentre shift for parabolic-like trajectories (p.e. 1%). Never- 
theless, this potential is not well-suited for studying more general 
cases (for instance it does not give the correct location of either r ms 
or r mb ). 



3.6 Accretion disc 

In this section we consider the simple model of the stationary, geo- 
metrically thin, and optically thick accretion disc first investigated 
by Shakura & Sunyaev ( 1973 1 in a Newtonian context and its exten- 
sion to a Schwarzschild spacetime by Nov ikov & Thorne] ( |1973| >. 
Under the present assumptions it turns out that the energy flux per 
unit area from the disc surface is independent of viscosity and is 
given by (see e.g. Fra nk et al.|2002| > 



D 



M 
4 7r r 



dr 



[h(r) - h(r in )} , 



(3.26) 



where M is the constant accretion rate and r in is the inner bound- 
ary of the disc at which the viscous stresses are assumed to vanish. 
In the standard thin disc model one takes r in = r ms . The total lu- 
minosity from the two faces of the disc is obtained after integrating 
Eq. J3T261 over the whole disc, i.e. 



L = 2 



D(r) 2-Krdr = eMc 



(3.27) 
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Figure 6. Energy flux per unit area D as a function of radius emitted from 
the surface of a geometrically thin, optically thick accretion disc. 



where e is called the efficiency of the accretion process. In the fol- 
lowing list we give the values of e for the different potentials: 

1 

12 



e N = — -E N (r ta )/c = 
e PW = —E rV j(ri„)/c z = 

E G (r in )/c 2 



: 0.083, 



e G 



1 

16 
1 

18 



0.062, 



0.056, 



(3.28) 



1 

18 



0.056, 



e s 



2V2 
~3~ 



0.057, 



from where we see that $ N w and $ G give an equally good ap- 
proximation to the relativistic value. Nonetheless, for the actual ra- 
dial dependence of D, the latter provides a better approximation 
(p.e. < 7.2%) than the former (p.e. 11%). In Figure|6] we com- 
pare D as obtained from general relativity with the results from the 
various potentials. 

3.7 Accretion inflow 



(3.29) 



Here we consider the analytic accretion model of Tejed a et al.| 
( |2012|[2013^ . In this model a rotating gas cloud of non-interacting 
particles accretes steadily onto a Schwarzschild black hole. In Fig- 
ure [7] we compare the streamlines of the relativistic solution with 
the ones obtained from the potentials $ N , <3? G , $pw, and <& N w for 
the following set of boundary conditions: 

M = 4M , 

M = 0.01M Q /s, 

r = 100 r g , 

f = -0.141c, 

r <j> = 0.038 c, 

0o = O. 

These boundary conditions are motivated by collapsing stellar 
cores leading to long gamma-ray bursts and were used in Tejeda 
|et al.||2012| l to make a comparison with one of the simulations of 
collapsar progenitors in Lee & Ramirez-Ruiz ( 2006 1. 

Note that the streamlines obtained from $ G are basically in- 
distinguishable from the exact relativistic ones. This is so because 
the energies of the incoming fluid elements are close to parabolic, 
i.e. £ ~ c 2 (for which case Eqs. (2T2| and (2T3} coincide with the 
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Figure 7. Streamlines of the steady-state accretion flow of a rotating gas cloud of non-interacting particles onto a Schwarzschild black hole. The boundary 
conditions are as given in Eq. l |3.29) . The figure shows a zoom-in of the first quadrant of the R-z plane (with R = r sin 8 and z = r cos 9). The black hole 
horizon is indicated with the dashed-line quarter-circle. 



exact relativistic equations). With a different choice of boundary 
conditions, the agreement with the relativistic results may not be as 
good as in Figure [7] although $ G is still in better agreement with 
the relativistic solution than $ PW , <&nw, or $ N . 

3.8 Tidal disruption 

In this section we apply the generalized potential to the tidal dis- 
ruption of a main sequence star by a supermassive black hole. To 
this end we have implemented the acceleration given in Eq. (A3) 
into the Newtonian SPH code that we had used previously to simu- 
late the tidal disruption of white dwarf stars by moderately massive 
black holes (Rosswog et al. 20091 and which has been described 
in detail in Ross wog et al.|p008} . For recent reviews of the SPH 
method see, e.g.[Monaghan (2005); Rosswog|( |2009[ >. 

Here we choose parameters that are identical to the ones of a 
general relativistic simulation run presented in Laguna et al. (1993 1: 
the masses are M B h = 10 6 M© for the black hole and M* = 
1 Mq for the star, the latter has a radius of R, = 1 R0 and is 
modelled by a polytropic equation of state with exponent T = 5/3. 
The star approaches the central black hole following a parabolic 
trajectory with an encounter strength /3 = r\/r p = 5, where r t — 
(M BU /M t ) 1/3 R* ~ 47 r E is the tidal radius and r p ~ 9.4 r E is 
the pericentre distance. As in the previous test cases we compare 
the results obtained with the generalized potential with those using 
"^n, $pw and $nw- 

The results for the different approaches are displayed in Fig- 
ure[8] For ease of comparison with La guna et al-H1993) , we simply 
show the trajectory followed by the centre of mass in each case to- 
gether with the particle positions projected onto the orbital plane at 
three different points of the trajectories. Obviously, for such a deep 
encounter relativistic effects lead to substantial deviations from the 
Newtonian parabola. In this case $ PW produces a qualitatively cor- 



rect result, though only about 70% of the relativistic pericentre ad- 
vance. As expected based on Fig.|5] $nw yields a pericentre shift 
with the opposite sign of general relativity, and therefore leads to 
a wider orbit than the purely Newtonian potential. Although the 
star passes the black hole with about 0.4 c, the hydrodynamic res- 
ult obtained from Q G is in very close agreement with the geodesic 
in Schwarzschild spacetime and the matter distribution closely re- 
sembles the one shown in Laguna et aL|( |1993fr (see their Figure 1, 
second column)F] 



4 SUMMARY 

We have derived a generalized Newtonian potential from the 
geodesic motion of a test particle in Schwarzschild spacetime in 
the limit where the velocity of the particle is much smaller than 
the speed of light. In spite of this, the present description provides 
an accurate approximation to the relativistic results even in cases 
where this condition is not entirely fulfilled. In addition to the 
standard Newtonian term, the generalized potential $ G includes a 
term that depends on the square of the velocity of the test particle 
and that can be interpreted as an additional gravitational attrac- 
tion from the central mass acting on the kinetic energy of the test 
particle. 

The new potential reproduces exactly several relativistic fea- 
tures of the motion of test particles in Schwarzschild spacetime 
such as: the location of the photon, marginally bound and margin- 
ally stable circular orbits; the radial dependence of the energy and 
angular momentum of circular orbits; the ratio between the orbital 
and epicyclic frequencies; the time evolution of parabolic-like tra- 
jectories; the spatial projection of general trajectories as function 

5 Bear in mind that for the derivation of <E>g we assumed that v <g c. 
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x/r g 

Figure 8. Tidal disruption of a solar-type star (M„ = 1 Mg , Ft* = 1 R©) by a supermassive black hole (M BH = 10 6 Mg). The figure shows the outcome 
of four runs of the same SPH code for each of the potentials $ N , &g, ^pw. and "3?nw- A common initial configuration is used in which the star is set at 
r = 200 r g along a parabolic trajectory with an encounter strength of /3 = 5. Only 10 4 SPH particles were used for each simulation. The figure shows the 
trajectory followed by the centre of mass of the star in each case (as projected onto the orbital plane), together with the projection of the SPH particles at three 
different points of the trajectory. For reference, we also show with a solid blue line the exact geodesic trajectory of a test particle in Schwarzschild spacetime. 



of the constants of motion and their pericentre advance. Moreover, 
the equations of motion derived from this potential reproduce the 
reversal of the centrifugal force at the photon orbit (see e.g. |Ab-| 
|ramowicz|19 90l. We are not aware that this relativistic feature has 
ever been reproduced before by any pseudo-Newtonian potential. 
Additionally, the obtained equations also provide a good approxim- 
ation for the time evolution of particles in free-fall, for the orbital 
angular velocity of circular orbits, and for the epicyclic frequencies 
associated to small perturbations away from circular motion (all 
these corresponding to the description made by observers situated 
far away from the central black hole). We have also applied $ G to 
the study of simple accretion scenarios, first for the thin disc model 
of Novikov & Thorne ( 1973 ) and then for the accretion infall model 
of Tejeda et al. ( 20 1 2 1, finding in both cases a good agreement with 
the exact relativistic solutions. 

As a further astrophysical application and a demonstration of 
the minimal effort required to implement the suggested generalized 
potential into an existing Newtonian hydrodynamics code, we have 
applied it to the tidal disruption of a main sequence star by a su- 
permassive black hole. For this, we implemented the equations of 
motion derived from $ G into the 3D SPH code described in |Ross-| 
|wog et aL"1 l|2008 ). The obtained results are in very good agreement 
with the relativistic simulation presented in Laguna et al. ( 1993 ). 

In Table [T] we summarise the accuracy with which various 
relativistic properties are reproduced by $ N , $ PW , $nw and $g. 



Potential 










At(r oc = 0) 


^ 16.7% 


s£ 23.8% 


< 9.3% 


exact 


At(r 00 = -0.3 c) 


^ 15.9% 


s£ 24.3% 


< 15% 


< 2.9% 


r,,h 




~ 33.3% 




exact 


r m b 




exact 


~ 13.4% 


exact 


r ma 




exact 


exact 


exact 


E c 


«: 50% 


< 12.5% 


< 3.3% 


exact 


h c 


^ 29.3% 


< 6.1% 


sS 29.3% 


exact 


n 


exact 


< 50% 


< 13.4% 


< 5.7% 


nil 


oo 


< 83.7% 


< 41.4% 


< 5.7% 


T 


^ 18.4% 


< 83.7% 


< 32.1% 


< 63.7% 


n - 7T 




~ 30% 


> 80% 


exact 


e 


45.7% 


9.3% 


2.9% 


2.9% 


D 


446.7% 


sc 58.9% 


< 11.1% 


< 7.3% 



Table 1. Maximum percent error obtained for the approximation of several 
quantities associated to the motion of test particles in Schwarzschild space- 
time with the four potentials <J> N , $ PW , $ NW and <J> G . In all the cases we 
have considered only r > 6r g . 



With the exception of the radial thrust T needed to keep a rocket 
hovering at a static position, we found that $ G provides a more ac- 
curate description of the motion of test particles in Schwarzschild 
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spacetime than any of the other potentials. For this reason, and 
given that a proper modelling of a realistic accretion scenario re- 
quires that (at least) all of these relativistic features are accurately 
reproduced, we consider the new potential as a simple, yet very 
powerful tool for studying accretion processes around a Schwarz- 
schild black hole. 
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APPENDIX A: ACCELERATION IN CARTESIAN COORDINATES 

The Cartesian coordinates (a;, y, z) are connected to the spherical ones [r, 6, (f>) in the usual way, i.e. 



x — r sin y cos < 
y = r sin 6 sin q 
z — r cos 9, 



r = \/ x 2 + y 2 + z 2 



tan 



1 



(j) = tan 1 (y/x) . 



from where we get 



xx + yy + zz = S ' a;* a;", 



r 4 ((9 2 + sin 2 # (j) 2 ) = (xy — y x) 2 + (x z ~ z x) 2 + (zy — y z) 2 — ( ^^e;jfc x 3 x 

i \ jk 

where x l — x,y,z and e ijk is the Levi-Civita symbol. Substituting Eqs. \P>..2) and \A3) into the Lagrangian in Eq. j2.1 1\ gives 



E«s < *'V + .£«(£i* c ** a!,i * 



r - 2r, 



r(r - 2r g ) 



+ 



GM 



from where we get the following expression for the acceleration components 

i.l/.r / 'J 1 ~ 



GMi" / 2r s 



2 r„ x z v— * 



r 2 (r — 2r g ) 



e jk ,x k x l 



(A.l) 

(A.2) 
(A.3) 

(A.4) 
(A.5) 



APPENDIX B: MISCELLANEOUS EXPRESSIONS 

In the following table we collect the different formulae associated to circular motion and plotted in Figures [2] and [3] The corresponding 
Newtonian expressions are also included to facilitate further comparison. 





Newton 


Schwarzschild 


Generalized Newtonian 


Paczyhski-Wiita 


Nowak-Wagoner 


T 

(h c ) 2 


— G y 


GM 


GM(r-2r E ) ! 


GM 


GA/(r 2 -6r g r+36 r 2 ) 


r 2 
GMr 


1 /r3(r-2r s ) 
GMr 2 


r 4 
GMr 2 


(r-2r g ) 2 
GMr 3 


r 4 

GM(r 2 -6r g r+36r 2 ) 


r-3r g 


r-3r g 


(r-2r g ) 2 


r 


E c 


GM 
2r 


GM(r-4r s ) 
2 r(r — 3 r g ) 


GA/(r-4r g ) 
2 r(r — 3 r g ) 


GM(t — 4r g ) 
2(r-2r g ) 2 


GM(r 2 -12r 2 ) 


Q 2 


GM 


GM 


GM(r-2r g ) 2 


GM 


GM(r 2 -6r g r+36r 2 ) 


r 3 




r 4 (r — 3 r g ) 


r(r — 2 r g ) 2 




(fi") 2 


GM 


GM(r-6r s ) 


GM(r-6r s )(r-2r g ) 2 


GM(r-6r fl ) 


GM(r 2 -36r 2 ) 


r 3 


r 4 


r 5 (r — 3r g ) 


r(r — 2 r g ) 3 


r° 
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